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Resolved images of self-gravitating circumstellar discs with ALMA 
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ABSTRACT 

In this paper we present simulated observations of massive self-gravitating circumstellar discs 
using the Atacama Large Millimetre/sub-millimetre Array (ALMA). Using a smoothed parti- 
cle hydrodynamics model of a 0.2 M disc orbiting a 1 M Q protostar, with a cooling model 
appropriate for discs at temperatures below ~ 160K and representative dust opacities, we have 
constructed maps of the expected emission at sub-mm wavelengths. We have then used the 
CASA ALMA simulator to generate simulated images and visibilities with various array con- 
figurations and observation frequencies, taking into account the expected thermal noise and 
atmospheric opacities. We find that at 345 GHz (870 jum) spiral structures at a resolution of 
a few AU should be readily detectable in approximately face-on discs out to distances of the 
Taurus-Auriga star-forming complex. 
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1 INTRODUCTION 

Circumstellar discs play an important role in the formation and evo- 
lution of both stars and planets, and as such have been the object 
of much study and observation in recent years. Millimetre wave- 
length surveys of star formi n g regions in Taurus dBeckwith et al.l 
ll99Ct iKitamura et alj 120021 : Undrews & Williams! 1 20051), Orion 
dEisner et alj|2008|) and p O phiucus dAndre & Montmerld [l994l : 
lAndrews & Willi ams 2007b) have provided extensive evidence of 
circumstellar discs, while optical images have been provided b y the 
Hubble Space Telesco pe (HST) of HH30 dBurrows et all 19961) and 
the Fomalhaut system l lKalas et"al1l2008h . However, relatively few 
systems have been imaged with great enough resolution to deter- 
mine the disc structure on scales of less than a few tens of AU. 

This may be about to change however, with the Atacama Large 
Millimetre/sub-millimetre Array (ALMA) due to come on line in 
the near future. With a minimum beam diameter of approximately 
5 milli-arcseconds at 900 GHz, ALMA should ideally provide res- 
olution down to » 2AU for discs observed in Orion, with sub-AU 
resolution for systems in Taurus-Auriga. It will therefore allow ob- 
servations of disc sub-structure in unprecedented detail. The dis- 
covery of structures within discs could have important implications 
for our understanding of the evolution of both the discs themselves 
dRice & Armitagell2009l) and the protostars they orbit. 

G aseou s circumstellar discs are expected to be turbulent 

dEbertlll994 lGammielll996l) . with the internal stresses that this 
turbulence induces being the principal driver of angular momen- 
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turn transport, and thus accretion. This turbulence will lead to sub- 
structure being present at all scales within the disc, with those on 
smaller scales (due to the magneto-rotational instability) remaining 
undetectable, while the larger scale structure induced by the gravi- 
tational instability should be resolvable with ALMA. 

It is expected that protostellar discs about Class 0/Class I 
objec t s will go through a self-gravitating phase dBertin & Lodatol 
1200 it IVorobvov & Basdl2005tlHartmannll2009l) as the infall rate 
from the gaseous envel ope is much greater tha n the accretion 
rate onto the protostar dVorobvov & Basul 1 20061) . In this phase 
large amplitu de spiral structures will form, driving accretion onto 
the protostar dLin & PringlJ 19871: lLaughlin & BodenheimeJll994l : 



lArmitage et alj|200ll : lLodato & Ricell2004 120051) and potentially 
leading to the formation of compan ion objects. These companions 
may vary in size from s t ellar m ass dBonnell|[T994l ; lBonnell & Bate! 
1 1994 IWMtworth eFafl 1 19951) to brown dwarf (BP) or plane- 
tary mass dRice et al .11 20051 : IStamatellos et al.|[2007l ; IClarkejl2009l : 
ICossins et al.ll2O10l) . It is possible that protoplanetary discs around 
Classic T-Tauri (Class II) stars may also ex hibit spiral patterns due 
to the presence of g ravitational instabi l ities dBertin &Lodatol2 001: 
iBolev et all 120061 : IVorobvov & Basj 120081 : ICossins et al .1 l2oToh . 
and indeed there are already possible detections of spiral str uctures 
in the discs of GSS 39 in Ophiuchus dAndrews et al.ll2009T) and in 
IRAS 16293-2422B dRodriguez et~ai1l2005l) . 

At radii greater than approximately 50AU, planet forma- 
tion through direct fragme ntation of these spiral overdensities into 
bound objects is possible dBossll 19971.1 19981: [ciai ke 2009; Rafikovl 
120091 : ICossins et alJl201Ch . While IWolf & D'Ang elo (2005]) have 
indicated that the forming giant (proto-)planets themselves may be 
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detectable using ALMA, the observability of the large-scale spiral 
structure itself within protostel lar and protoplanetary discs, though 
implied jTesti & L eurini 2 008h . has remained undemonstrated. In 
this paper we therefore present a simple self-gravitating disc simu- 
lation, and from it we derive mock observations of disc systems 
at the resolutions and sensitivities that should be possible with 
ALMA. Hence in section [2] we detail the simulation, and in sec- 
tion[3]we discuss how the mock observations are generated from it, 
taking into account telescope effects and sensitivities. In section [4] 
we present the observations for various system and telescope pa- 
rameters, and finally in sectionf5]we discuss the significance of our 
results. 



in our simulation will have qualitatively similar structure forma- 
tion, with the uncertainty lying primarily in the amplitudes of the 

density perturbations. 

Furthermore, in ICossins et af] J2010h we characterise the be- 
haviour of flr coo i with radius, and show that for accretion rates 
~ 1CT 7 M yr~' (typical for Classic T Tauri objects) at radii of 
20 - 50AU the value of Q.t caol is ~ 10 1 - 10 2 , decreasing with 
increasing radius. Hence although the cooling formalism given in 
equation [2] is very simple, it produces spiral structure in the cor- 
rect modes and at approximately the correct amplitudes for the 
outer radii of discs, as long as we choose f} in the range > 10. It 
is therefore useful as a means of generating test cases to investigate 
whether such structures would actually be observable, particularly 
in the outer parts of discs. 



2 NUMERICAL SIMULATION OF STRUCTURE 
FORMATION 

Gravitational instabilities in discs have been exte nsively studied 
numerically using both grid-based Eulerian codes |Gammiell20o"ll ; 



IPickett et al.ll2003l ; iBolev et alj|2006l: 



2004 120051: iNavakshin et al j |2007|; Alexander et al.l 120081: 



Kratter et al 



20101) and La- 



grang ia n part ic le-based methods dRice et al. 20031; lLodato & Rice] 



Forgan & Ri ce 2009 ; Cossins et al.l2009l) and the reader is referred 



to these for further details of the dynamics of gravitationally unsta- 
ble discs. The simulation we use to gen erate the mock observ ations 
was performed in the manner detailed in lCossins et al. using 
a 3D smoothed particle hydrodynamics (sph) code, a Lagrangian 
hydrodyn amics code capable of modelling self-gravity (see for 
exam ple, iBenj Il990l : iMonaghanl 1 19921 120051 : IPrice & Monaghanl 
1 20071) . Our system consists of a single point mass orbited by 
500,000 sph gas particles, with the central object free to move 
under the gravitational influence of the disc. All particles evolve 
according to individual ti me-steps governed by the Courant 
condi tion, a force con dition dMonaghar] 1 1 99% and an integrator 
limit dBateet all 19951) . 

We allow the discs to cool towards gravitational instability by 
implementing a simple cooling law of the form 



diij 
dt 



(1) 



where u t and f coo] ,- are the specific internal energy and cooling time 
associated with the i' 1 ' particle. We then determine the cooling time 
f coo i through the simple prescription that 



(2) 



where fl is the angular frequency and /? is a fixed parameter 
throughout each simulation. Although this is clearly an ad hoc cool- 
ing function, it can be used as a simple param eterisation in order to 
conduct controlled numerical experiments. In lCossins et al. 1 120091) 
we found that for a given disc mass, the spiral structures formed 
through the gravitational instability (as characterised by the radial 
and azimuthal wavenumbers of the excited modes) are independent 
of the cooling, but that the strength of the modes (characterised 
by the relative RMS amplitude of the surface density perturbations 
(iSE/E)) varies such that 



\5L\ 



(3) 



Therefore, notwithstanding the simplicity of our cooling prescrip- 
tion, we may reasonably assume that the spiral structures formed in 
physical systems with characteristic masses corresponding to those 



2.1 Simulation Details 

The simulation used to generate the mock observations described 
here consists of a 0.2 M disc about a 1.0 M Q mass star, ex- 
tending out to approximately 25 AU. Altho ugh this is a relatively 
high mass ra tio, it is within observed bounds l lAndrews & Williams! 
l2005Ll2007bl) . and is plausible for the early (Class I) stages of inter- 
mediate mass pr otostellar evolution, when the disc is expected to be 
self-gravitating dBertin & Lodato! 1200 ll ; IVorobvov & Basul 120051: 
lHartmaniJl2009l) . 

Our initial conditions consist of a disc of gas particles on cir- 
cular orbits, distributed with a surface density profile S ~ R~ y 
with y = 1.5, as for the Minimum Mass Solar Nebula (MMSN, 
IWeidenschillingl f"l977l) . Note that observational constraints for y 
range f rom 0.4 < y _g L0 for protoplanetary discs in Ophi- 



uchus 1 Andrews et alj 20091) to 0.1 < y < 1.7 in Taurus 



l lAndrews & W illiam sl2007al) . so this value is within observed lim 
its. 

The disc is initially in approximate vertical hydrostatic equi- 
librium with a Gaussian distribution of particles where the scale 
height H = c s /f2 and c s is the sound speed. The azimuthal ve loci- 
ties take into account both a pressure correction dLodatol2007l) and 
the enclosed disc mass and any variation from dynamical equilib- 
rium is washed out on the dynamical timescale. The initial temper- 
ature profile is such that c 2 ~ R~' 12 , with the minimum value of the 
Toomre parameter Q mm = 2 occurring at the outer edge of the disc. 
In this manner the disc is initially gravitationally stable throughout 
its radial range. 

Note that the sph code and the initial cond i tions used ar e ex- 
actly the same as were us ed in lLodato & Ricel d2004l 12051) and 
ICossins et alj d2009l |2010|) . excepting the fact that here we have 
used a mass ratio of 0.2, and the reader is directed to these for fur- 
ther details. Finally, in terms of cooling, we set j} = 7 and use the 
cooling formalism described in equation [T] above. This is approxi- 
mately in the expected range, and is low enough to avoid spurious 
nume rical heating effects due to artificial viscosity dLodato & Rice! 
l2004h . 



2.2 Disc Evolution 

Although the thermal profile of the disc initially ensures it is gravi- 
tationally stable at all radii, it is however not in thermal equilibrium. 
As the simulation evolves, the disc therefore cools towards gravi- 
tational instability, which is initiated when Q ~ 1, after approxi- 
mately 1000 years. The disc then settles into a marginally stable, 
quasi-steady dynamic thermal equilibrium state, characterised by 
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Figure 1. Simulated surface density perturbations in a 0.2 M Q disc about a 
1 .0 M Q protostar. The gravitationally induced spiral waves that impart heat 
to the disc are clearly visible. 



the presence of spiral density waves that propagate across the face 
of the disc, and which provide heat (through shocks) to balance the 
imposed cooling. These spiral waves are clearly seen in the surface 
density, as shown in Fig. [T] The simulation was run for approxi- 
mately 10,000 years (equivalent to ~ 10 cooling times at the disc 
outer edge), suggesting that during the self-gravitating period of 
protostellar evolution (expected to last a few x 10 5 years during and 
immediately after the infall phase) circumstellar discs should in- 
deed be able to reach this quasi-steady state. Since the surface den- 
sity varies on the longer viscous timescale, the original X ~ R~ 3/1 
profile remains largely unchanged throughout the simulation. 

Once the disc reaches the self-regulated dynamic thermal 
equilibrium state, the disc temperature falls to approximately 20 - 
40 K, with the minimum at the disc outer edge, as shown in Fig. [2] 
Between roughly 10 and 25 AU the temper ature falls off as R [ , 
in reas onable agreement with observations ^Andrews & Williams! 
l2007al) . Note that the slight temperature rise outside of R « 25 AU 
is due to very low density, higher temperature gas external to the 
main bulk of the disc. 



3 GENERATION OF MOCK OBSERVATIONS 

Having run the hydrodynamic simulations, it is then necessary to 
use the disc parameters to create flux maps of the objects as they 
may appear using ALMA. Throughout the following we assume 
that the disc is face-on to the observer. 

In order to calculate the emission we use the individual par- 
ticle densities pi generated from the sph simulations, and calculate 
the particles' absorption coefficients a vi at frequency v using 



(4) 



where k v is the opacity of the disc at frequency v. 

For a face-on disc the optical depth t v at frequency v is defined 
in the following manner 



£ 



a y (z) dz, 



(5) 



and hence for a disc of SPH particles, we may evaluate the ver- 
tical optical depth at any point on the disc face by approximating 
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Figure 2. Azimuthally averaged temperature profile of the disc once it has 
settled into the dynamic thermal equilibrium state. The dashed line shows 
an profile. 



this integral as a sum over all the relevant interacting particles. We 
therefore use the following approximation for the optical depth of 
the disc as seen by a distant observer 



/ \ PiKylt 



(6) 



where w\ is a weighting function related t o the particl e's mass, den- 
sity, and the sph smoothing kernel (see IPricd [2007h and where i 
denotes all particles along a given line of sight through the disc. 

We assume that radiation from the disc is in thermal equilib- 
rium with itself and thus that the disc emits as a black body, with 
the source function S v at frequency v given by the Planck function, 

S V = B V (T) = ^f- } (7) 
c 2 e hv > kT - 1 

where h is Planck's constant, c is the speed of light in vacuo, k is 
Boltzmann's constant and T is the temperature of the emitter. For 
simplicity, we assume that our discs are vertically isothermal with 
temperature T, (which we obtain via a vertical average from the 
simulations) meaning that at each point in the disc the source func- 
tion is constant with height, and is given by B V (T). Modelling the 
disc in this manner, as a geometrically thin structure with no verti- 
cal temperature gradient, is obviously a very crude approximation. 
Nevertheless, since we are only interested in predicting the emis- 
sion in the (sub-)millimetre range, this approximation is justified 
as the bulk of the emission will indeed come from the relatively 
isothermal layer located about the disc midplane. 

For such a constant source function, the specific intensity (or 
surface brightness) l v at frequency v and optical depth r v is given 
by 



I v = B v (T)(l-e- T >). 



(8) 



Note that this formulation for the surface brightness is often used 
in reverse to in fer disc properties from the (sub) millimetre s urface 
brightness, e.g. lBeckwith et all i ll 990h : I Andrews et al. I J2009t) . Fur- 
thermore, it is also clear that once the disc becomes optically thick 
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its emission is determined solely by its temperature, whereas in op- 
tically thin regions the emission is exponentially dependent on t v . 
Optically thick structures forming in optically thin parts of the disc 
are therefore likely to show greater variation in intensity (and thus 
be more readily observable) than structure in purely optically thick 
regions. 



3.1 Dust opacities 

A critical parameter in the above calculation is the opacity of the 
disc k v . At temperatures belo w x 160K the Rosseland mean opacit y 
becomes dominated by ices jBell & Linll99ilCossins et alj201Ch . 
but the specific value of the opacity at a given frequency k v is det er- 
mined by the grain size distribution jMivake & N akagawa 19931) . 

To determ i ne the value of k v we use the power-law model of 
iBeckwith et al. I fl99ch . such that 



1V12 



(9) 



where v l2 = 10 12 Hz, and where K i2 is the value of the opac- 
ity at this fiducial frequency. The normalisation constant K12 is 
relat ively poorly constrain ed, with values ranging from 0.1 cm 2 



g _1 (Beckwith et ajJ |l990h to approximately 0.016 l lKramer et al .1 
" 20091) 



1 1998 ; Rafikov 2009h . For protoplanetary and protostellar discs the 



power-law index p can be determined from the disc SED in the 
mm/sub-mm range, and is found to be in the region 0.3 < < 1.5 
dKitamura etakl I2OO2I ; iTesti et all 120031 : iRicci et alj 120091) . Note 
that this rang e is below that expected for the interstellar medium 
p » 1.7-2 dFinkbeiner et alj| 19991 : IChakrabarti & McKedl200& 
lHartmannll2009l) . which is attributed to grain processing within the 
disc. 

For our models we use a fiducia l dust o pacity k 12 = 0.025 cm 2 
g _1 , and in accordance with Ra fikovl J2009h we assume throughout 
that/? = 1.0, such that 



0.025 



(l0 12 Hz) 



cm- g 



(10) 



3.2 The ALMA Simulator 

To produce realistic simulated ALMA observations of the models 
we have used the ALMA simulator simdata in CASA0. The ver- 
sion of the simulator we used (2.4) allows us to use the latest ver- 
sion of the planned ALMA antenna configurations, the expected re- 
ceiver noise based on technical specifications and the contribution 
due to the atmosphere, itself based on input values for the atmo- 
spheric temperature T atm and the optical depth at the frequency of 
the simulations. In all the simulations presented here we assume 
r atm = 265 K, and th e optical depth is computed using the ATM 
atmosphere models of lPardo et al.l J2002h . We use typical Chajnan- 
tor conditions and an amount of precipitable water vapour (pwv) as 
expected for dynamical scheduling of the observations. In Table [T] 
we therefore show the water content and opacity of the atmosphere, 
the expected theoretical noise and the angular and linear resolution 
of the simulations for discs at 50 pc and 140 pc for each frequency. 
Note that we have also run simulations for the lowest (45 GHz) 



1 CASA (Common Astronomy Software Applications) is being primar- 
ily developed by ALMA and the National Radio Astronomy Observatory 
(NRAO) as the main off-line data reduction package for ALMA and EVLA, 
the Expanded Very Large Array, http : //casa . n rao . edu| 



and highest (870 GHz) planned ALMA bands - contrary to the in- 
termediate frequencies, these frequency bands will not be available 
when ALMA first comes on-line. High frequencies receivers are 
planned to be introduced during the early years of ALMA science 
operations, while the low frequency band is under discussion for 
the ALMA development programme. 

The array configuration used varied for each frequency, as 
we aimed to use the configuration that offered the best compro- 
mise between surface brightness sensitivity and angular resolution. 
Note that at low frequencies the angular resolution is limited by 
the largest available array configuration. In all cases we performed 
aperture synthesis simulations with transit duration of 2 hrs. 



4 RESULTS 

Using the disc simulation shown in Fig.Q] we investigated the emis- 
sion at the various frequencies shown in Table [T] Based on the 
requirement that the optical depth should vary between optically 
thick in the spiral arms and optically thin in the inter-arm regions 
(to maximise the contrast in emission across these regions) we have 
principally considered the results at 345 GHz. Furthermore, at this 
frequency we expect the resolution to be roughly 1 AU at both TW 
Hydrae and Taurus-Auriga distances, while the telescope sensitiv- 
ity should not be a limiting factor. To illustrate the effects of varying 
the frequency, the optical depth at 45, 345 and 870 GHz is shown 
in Fig.[3]for comparison. Also clear from this figure is the expected 
increase in optical depth with frequency due to increasing opacity. 



4.1 Simulated ALMA Images 

Having generated the specific intensity map by the method de- 
scribed in section [3] we have used them to simulate ALMA ob- 
servations as described in section [3~2l In Figs[4]and[5]we therefore 
show the results of our simulations for discs at a distance of 50pc 
(TW Hya) and at 140pc (Chamaeleon, Ophiuchus, Taurus-Auriga). 
Note that our simulations only include the effects of thermal noise 
from receivers and the atmosphere, but do not take into account 
calibration uncertainties and residual phase noise after calibration. 
These effects are likely to be most important at high frequencies 
and long baselines, so the simulated maps at 680 and 870 GHz, es- 
pecially for the 140pc case, represent observations carried out in 
ideal conditions and with excellent calibrations. 

The simulations show that the predicted spiral structure is 
readily detectable at all simulated frequencies at the 50 pc dis- 
tance of the TW Hya association. In the case of star-forming re- 
gions at 140pc the situation is less clear cut. At low frequencies 
(< 100 GHz) even ALMA will probably not provide the angular 
resolution required to image the spiral structure clearly, whereas at 
the highest frequencies, as noted above, the simulations are prob- 
ably over optimistic. Nevertheless, our simulations show that at 
220 and 345 GHz (ALMA Bands 6 and 7), the predicted structures 
should be clearly detectable. 

Finally, in Fig. [6] we show the predicted observability of a disc 
at 410pc (equivalent to the Orion Nebula Cluster distance) imaged 
at 345 GHz. While the non-axisymmetric nature of the disc is clear, 
the spiral structure per se is not well resolved, and thus we infer 
that even with ALMA such structures will not be conclusively de- 
tectable. 
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Figure 3. Comparison of the (logarithmic) optical depth across the disc face at frequencies of 45 GHz (left), 345 GHz (centre) and 870 GHz (right). The 
greatest contrast between optically thick and optically thin regions across the arm/inter-arm regions in the bulk of the disc is to be seen in the 345 GHz case. 
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0.05 


3 


120 


6 


120 


16.8 


100 


2.3 


0.03 


4 


50 


2.5 


50 


7 


220 


2.3 


0.1 


10 


60 


3 


20 


2.8 


345 


1.2 


0.2 


20 


40 


2 


25 


3.5 


680 


0.5 


0.6 


60 


50 


2.5 


35 


4.9 


870 


0.5 


0.7 


100 


60 


3 


45 


6.3 



Table 1. Atmospheric conditions, expected sensitivities and angular and linear resolutions as a function of frequency for the simulated observations, assuming 
distances of 50 pc (corresponding to the distance of the TW Hya association) and 140 pc (roughly corresponding to the Taurus-Auriga, Ophiucus and 
Chamaeleon star forming regions). The amount of precipitable water vapour (pwv, in mm) has been chosen according to the current ALMA dynamical 
scheduling expectations. Note that at the higher frequencies the varying resolutions at different distances is due to the fact that we have considered different 
ALMA configurations in order to obtain an optimal compromise between resolution and sensitivity. On the other hand, at low frequency, the resolution is 
dictated by the largest available array configuration, and thus remains constant with distance. 
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Figure 4. Simulated aperture synthesis ALMA images for a disc at 50 pc, with a transit duration of 2 hours. From left to right and top we show simulations 
computed for an observing frequency of 45, 100, 220 (top) and 345, 670 and 870 GHz (bottom). Axis scales are in arcseconds. Contours start at 0.01 and are 
spaced by 0.08 mJy/beam at 45 GHz, start at 0.08 and are spaced by 0.2 mJy/beam at 100 GHz, start at 0.5 and are spaced by 0.5 mJy/beam at 220 GHz, start 
at 0.8 and are spaced by 0.8 mJy/beam at 345 GHz, start at 2 and are spaced by 2 mJy/beam at 680 GHz, start at 4 and are spaced by 4 mJy/beam at 870 GHz. 
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5 DISCUSSION 

In this paper, we have used a 3D, global SPH simulation of a mas- 
sive (0.2 M Q ) compact (R out = 25 AU) self-gravitating disc about 
a young star (1.0 M Q ) to demonstrate that the spiral modes excited 
by the gravitational instability should be detectable in face-on cir- 
cumstellar discs using ALMA. At distances comparable to the TW 
Hydrae association (~ 50 pc), such spiral density waves are read- 
ily apparent with observation times of 2 hours, whereas at Taurus 
distances of ~ 140 pc a careful choice of the observing frequency 
and excellent observing conditions may be required for significant 
detections. Our results suggest that structure in such discs in Orion 
(~ 400pc) will most likely not be resolvable. 

In order to generate these predicted observations we have used 
temperature and density maps provided by numerical simulations, 
together with an empirical relationship for the dust opacity of cir- 
cumstellar material to obtain the optical depth at each part of the 
disc face. In a precisely similar manner to that used to obtain disc 
properties from sub-mm observations dBeckwith et alj[l990h we 
have then been able to determine the specific intensity of the disc 
emission across the disc face, which has then been used as input for 
the ALMA simulator in CASA to produce simulated observations. 
We find that observation frequencies of around 345 GHz (870 /jiti) 
are ideal for this kind of observation, as we find the spiral arms to 
be optically thick, whereas the inter-arm regions are optically thin, 
maximising the contrast in emission between the regions. It should 
be noted however that this 'ideal' frequency range is dependent on 
our assumptions concerning the dust opacity, in which there is con- 
siderable scatter. 

There are certain limitations to our model which should be 
noted however. Our stellar and disc masses are both at the up- 
per end of the expected distr ibutions ^Andrews & Williams|[200a, 
l2007bl : iBeckwith et alj[l990h , although in the early stages of star 
formation (Class I objects) such high disc masses are not unrea- 
sonable due to the infalling envelope. Of necessity, in order to be 



self-gravitating our discs are cold (20 - 40 K), which implies both 
a relatively low ambient temperature (not unreasonabl e, since giant 
molecular clouds tend to have temperatures of ~ lQjC. lMvers et al.l 
ll983l ; lMvers & Bensonll 19831) and for heating from the protostar to 
be negligible. While this latter assumption is clearly unlikely to be 
valid for the surface layers of the disc that are irradiated directly by 
the star, the disc midplane (which dominates the emission) is likely 
to be cold enough to justify this assumption ^Andrews & Williams! 
120051 : iDullemond eTal]|2007l ; iD'Alessio et al.ll 19981) . In a similar 
manner, this assumption of a colder inner layer allows us to ign ore 
the effects of the magneto-rotational instability dGammiell 19961) . 

Given that our disc is quite compact, and many discs are 
obs erved to extend out to much larger radii (~ 10 2 — 10 3 
AU. I Andrews et alj|2009l: lEisner et al.ll2008l: I Andrews & Williams] 
l2007bl : lKitamura et al.l2002l) the colder regions outside of ~ 50 AU 
are if anything more likely to show evidence of gravitational insta- 
bility that at the radii we simulate, and indeed any spiral structures 
forming at large radii will be more easily resolvable. In this sense 
therefore our predictions may even be conservative in terms of the 
maximum distance at which spiral structures may be detectable. 

Likewise our cooling prescription, although simplistic, is valid 
for regions at large radii where the temperature is below the ice sub- 
limation temperature, and is therefore not an unreasonable simpli- 
fication. It should be noted that the ratio of the cooling time to the 
dynamical time /? = t coa] Q determines the amplitude of the spiral 
perturbation and hence the contrast in the simulated images. The 
value p = 7 adopted in this paper is in the right range for discs at a 
few tens of AU. However, a larger value of /? (that is, a less efficient 
cooling) would provide a relatively smaller contrast in the ALMA 
images. 

Finally, note that in this paper we have considered the contri- 
bution to the sub-millimetre emission due solely by the disc. In the 
earliest phases of star formation, the system might show a substan- 
tial emission on larger scales, produced by the infalling envelope 
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feeding the disc. This larger scale contribution has been neglected 
in the present paper. 

As noted in the Introduction there are already possible de- 
tections of spiral stru ctures in the discs of GSS 39 in Ophiuchus 
dAndrews et alj|2009l) and in IRAS 16293-2422B (Rodriguez et al. 
120051) . However the structures in GSS 39 are not robust at the 
3cr level dAndrews et aHl2009l) . and those in IRAS 16293-2422B, 
whilst appearing to be genuine, may plausibly be due to inter- 
actions with a companion. Furthermore, the presence of plane- 
tary mass objects within the disc may als o lead to the excitation 
of sp i ral density waves (see f or instance Goldreich & Tremaind 
Il98d : lQgilvie&Lubowll2002l : lArmitage & Ricdl2005h . However 
in this latter case the one-armed spiral wa kes are in general tightly 
wound (dependent on the disc thickness. lOgilvie & Lubow 2002) 
and would therefore be difficult to detect directly. Furthermore, 
analysis by Wolf & D'Angelo] d2005l) on the observability of em- 
bedded gas giant planets using ALMA suggests that the spirals 
would only be visible very close to the planet itself. Any detec- 
tions of large scale spirals in isolated protostellar discs are there- 
fore likely to be due principally to the action of the gravitational 
instability. 

Confirmed, unambiguous observations of gravitationally in- 
duced spiral structures within protostellar discs would be valu- 
able for a number of reasons. As the gravitational instability is 
expected to operate during the early phases of star formation, pro- 
cessing the infalling envelope and allowing rapid accretion onto 
the protostar, such detections would validate this mechanism for 
growing the masses of protostars. More controversially, it may en- 
able models of brown dwarf/lo w mass stellar companion forma- 
tion th rough disc fragmentation 1 Bate et al J 20021: Stamatellos et al.l 



120071 : IStamatellos & Whitworth] l2009l : IClarkei (20091) to be vali- 
dated, as the presence (or otherwise) and amplitudes of spiral arms 
would allow constraints to be placed on the numbers of such com- 
panions that may be expected. 

In a similar vein, detections of spiral features may enable us 
to determine the dominant mode of planet formation at various 
radii, about which there is some debate. While the s t andard core- 
accretion model is favoured at low radii dBolevll2009l ; lKlahrll2008l ; 
iBodenheimer et al]|2000l : lLissauedll993l) . direct fragmentation of 
gravitationally induc ed spirals remains a candidate mechanism at 
radii above - 50 AU dBosslll997l ; lBolev et alJl2006l ; lRafikovll2009l ; 
ICossins et al.ll201ol) . As the presence of a peak in the 1.3 cm emis- 
sion around HL Tau has been put forward as a promising candidate 
for a planetary mass companion formed through gravitational insta- 
bility, detections of the spiral wave progenitors of such companions 
would provide significant backing for this mechanism. Moreover, 
the presence of large amplitude spiral density perturbations may be 
important for the formation and growth of planetesimal s, both due 
to the concentration of the dust fract ion within the arms l lRice et all 
|2004 120061 : IClarke & Lodatoll2009l) and further due to the possi- 
ble sc attering of planetesimals by the spiral potential dBritsch et al.l 
120081) . In either case, observations of the spiral arms themselves 
would place constraints on, and therefore allow us to discriminate 
between, the two planet formation modes at large radii. 
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Figure 6. Simulated 345 GHz aperture synthesis image of a disc at 410 pc, 
with a transit duration of 2 hours. Contours start at and are spaced by 60 
mJy/beam. Clear asymmetries are present, but the underlying spiral struc- 
ture is not well resolved. 
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